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Abstract 

We study the stationary and axisymmetric non-convective differentially rotat- 
ing perfect-ffuid solutions of Einstein's field equations admitting one conformal 
symmetry. We analyse the two inequivalent Lie algebras not exhaustively consid- 
ered in [|J and show that the general solution for each Lie algebra depends on one 
arbitrary function of one of the coordinates while a set of three ordinary differ- 
ential equations for four unknowns remains to be solved. The conformal Killing 
vector of these solutions is necessarily homothetic. We summarize in a table all 
the possible solutions for all the allowed Lie algebras and also add a corrigendum 
to an erroneous statement in concerning the differentially rotating character 
of one of the solutions presented. 



*Also at Laboratori de Ffsica Matematica, IEC, Barcelona. 
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1 Introduction 



In a recent paper QXJ] we considered the study of stationary (non-static) and axisymmet- 
ric non-convective differentially rotating perfect-fluid solutions admitting one proper 
conformal Killing vector. A purely geometric result [0 for stationary and axisymmet- 
ric spacetimes (with a well-defined axis of symmetry) admitting one conformal motion 
implies that both the stationary and the axial Killing vectors necessarily commute with 
the conformal Killing. Thus, the Lie algebra of the three dimensional conformal group 
can only adopt four inequivalent forms, which were labeled in |T[ as Abelian Case and 
Cases I, II and III, and which correspond to 



Abelian Case 


Ik 


= 0, 


Case I 


Ik 


= bk, 


Case II 


Ik 


= <£, 


Case III 


Ik 


= Vf}, 



— » — * 

where £ is the timelike Killing vector, fj is the axial Killing vector, k is the conformal 
Killing and b, c and v are arbitrary non-vanishing constants (these constants can be 
normalised to one, but as they usually carry dimensions we prefer to retain them). In |]J 
the general solution of Einstein's field equations for the Abelian and Case I was found. 
The general solution for the Abelian Case depends on an arbitrary function of a single 
variable and the perfect fluid satisfies the barotropic equation of state p = p+const. The 
solution is Petrov type D and the fluid velocity vector lies in the two plane generated 
at each point by the two repeated principal directions of the Weyl tensor. This solution 
was first found under completely different hypotheses by one of us In Case I, the 
general solution is an explicit Petrov type D metric with the fluid velocity vector lying 
outside the two-plane generated by the repeated principal directions and the perfect 
fluid satisfying the barotropic equation of state p + 3p = 0. Let us remark here that 
|l[ contained an erroneous statement concerning this solution which is corrected at 
the end of this paper: the fluid velocity vector for this solution was claimed to rotate 
differentially when the solution is, actually, rigidly rotating. Then, as a consequence of 
the results in [H , this solution is contained in the general Wahlquist family |J of Petrov 
type D rigidly rotating perfect fluids with equation of state p + 3p =const. 

The Cases II and III were also considered, but the study of solutions was only 
performed under the severe restriction imposed by an Ansatz of separation of variables 
in one of the metric functions. For the first two cases in ([!]) this separation of variables 
is a consequence of one of the Einstein field equations but, for the remaining Cases II 
and III, it had to be imposed as an additional assumption in order to handle the field 
equations. The solution under this assumption was found and the conformal Killing 
vector turned out to be, in fact, homothetic. 
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The aim of this paper is to complete the study performed in ]I[ by analysing the 
general non-convective differentially rotating perfect-fluid solution of Einstein's field 
equations for a stationary and axisymmetric spacetime with one conformal symmetry 
such that the Lie algebra of the conformal group belongs to one of the Cases II or III 
above. It turns out that there exists one perfect-fluid family of solutions for both cases. 
The dependence of all metric coefficients of these solutions on one of the coordinates 
is explicit. The dependence on the other non-trivial coordinate can be obtained by 
solving a remaining compatible system of ordinary differential equations. The general 
solution of this system depends on an arbitrary function which cannot be reabsorbed 
by any coordinate change. It also results that the conformal Killing for these two 
families is necessarily homothetic which proves that no non-convective stationary and 
axisymmetric perfect fluid admitting one proper conformal motion exists apart from 
those in We remark that the families we will find in this paper contain as 

particular cases the solutions for Cases II and III considered in |TJ. 

The plan of the paper is as follows. In section 2 we recall the line-element and the 
form of the conformal Killing for the Lie algebras II and III. We rewrite the metric in 
a way which allows us to consider the two cases together. Then we write the Einstein 
field equations for a perfect-fluid energy-momentum tensor and concentrate our efforts 
in a subsystem constituted by two partial differential equations for a single function 
ty(x, y) (the conformal factor of the metric). In section 3 we find the general solution 
of this system of equations which turns out to be a set of four different families that we 
label as Solutions 1, 2, 3 and 4. Section 4 is devoted to impose the remaining Einstein 
equation (which is quadratic in the second derivatives) for each of the four solutions. 
We prove that Solution 1 does not represent a perfect fluid. Regarding Solutions 2 and 
3, we show that they can be treated together and that the quadratic equation gives rise 
to a system of ordinary differential equations which can be shown to give no solution 
with a regular axis of symmetry. In an Appendix we proof this statement completely 
for the Case III and give the main indications regarding the Case II (the proof for this 
case is much longer but still not difficult to perform). Let us remark here that the 
assumption of regular axis is not an additional one in our work because the restriction 
to the four possible Lie algebras (^) is based essentially on this regularity of the axis 

Then we consider Solution 4, which makes the quadratic equation compatible thus 
giving perfect-fluid solutions of Einstein's equations which are analysed in Section 5, 
where we rewrite the line-elements for these solutions in the simplest forms and write 
down the remaining ordinary differential equations including a discussion on the com- 
patibility of the system. We also give the expressions for the density and pressure 
and for the fluid velocity vectors. We finish the section by enclosing a table which 
summarizes all the non-convective stationary and axisymmetric perfect-fluid solutions 
of Einstein's equations admitting one conformal Killing vector. We end this paper by 
including a corrigendum for an erroneous statement in |TJ concerning the solution in 



3 



Case I (see above). 



2 Canonical forms of the metric and Einstein field 
equations 

As we have stated in the Introduction, we are interested in studying the Einstein field 
equations for a non-convective, differentially rotating perfect fluid in stationary and 
axisymmetric spacetimes admitting one conformal motion such that the conformal al- 
gebra belongs to Cases II or III defined above. The canonical form of the metric in 
these two cases is explicitly given in ffl] where it is shown that there exist coordinates 
in which the line-elements and the conformal Killing vectors take the forms 



Case II 

ds 2 



^ 2 (x,y) 



-F(x) e~ cy dt + P(x)d(f) + 



Q\x) 



F{x) 



+ dx 2 + dy 2 



(2) 



k-ct— + — 
dt dy 



Case III 



ds 2 



^ 2 (x,y) 



-F(x)dt 2 + 



F(x) 



+ (P(x) - vy) dt + dx 2 + dy 



(3) 



k - vt— — 

d(j) dy' 

In both metrics the dependence of the function ty(x,y) on the variable y must be non- 
trivial because, otherwise, the conformal Killing k becomes a true Killing field. 

As we will see in what follows, it is very convenient to study the two Cases II and 
III together as they have similar properties. It is obvious that the two line-elements 
can be fused into the following expression 



ds 4 



^ 2 (x,y) 



-F(x)e- 2cy dt 2 + ( d<P + e~ cy (P(x) - vy) dtf + dx 2 + dy 2 

F(x) v ' ' 



(4) 



where the expressions relating the functions F(x), P(x) and Q{x) with F(x), P{x) and 
Q(x) for Case II can be easily obtained. Case III is recovered by setting c = and Case 
II corresponds to v — 0. Thus, we can set the product cv = wherever it occurs in the 
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calculation of the Einstein tensor. The orthonormal tetrad that we will use in order to 
write the components of the tensor quantities is given by 

0° = -VFe- cy dt, 1 = --Q=( d(f> + e~ cy (P - vy) dt] , 2 = -dx, s = —dy. 

q> q> y/p v / q? q? 

We are interested in perfect-fluid solutions of the Einstein field equations, so that the 
energy-momentum tensor takes the following standard form 

T a/ 3 = {p + p) u a U/3 + pg a(3 , 

where p and p are the energy density and the pressure of the fluid respectively and u 
is the fluid velocity vector. The fluid motion is taken as non-convective, that is to say, 
the velocity vector u lies on the two-plane generated at each point by the two Killing 
vectors £ and fj. In the above orthonormal tetrad the fluid one-form u is given by 

u = uq6° + uiO 1 , lip — u\ = 1, 

so that the Einstein field equations, in c = 8irG = 1 units, read 

S 00 = (p + p)u 2 -p, S i = (p + p) u ui, S n = (p + p)ul+p, 

S22 = 5*33 = p, Sq2 = 5*03 — $12 — 5*13 = S23 = 0, 

where S a p stands for the components of the Einstein tensor in the {0 a } cobasis. Some 
of these equations can be rewritten involving only the Einstein tensor, while the energy 
density, pressure and velocity vector are calculated once these equations are solved. Due 
to the separation of the metric into two orthogonal blocks, some of the Einstein tensor 
components are identically vanishing in the tetrad we are considering. Specifically, it 
can be shown by direct calculation that 

Sq2 = 0, 5*03 = 0, S12 = 0, S13 = 

and consequently, the only non-trivially satisfied Einstein equations are 

(•S'oo + S22) (Sn — S22) — S 01 = 0, (5) 
S 22 - S33 = 0, (6) 
S23 = 0. (7) 

We will now concentrate on the last two equations which read (after dropping 

some global non- vanishing factors) 

d xy q> = H(x)% d xx q> - d yy q> = G(x)q>. (8) 
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where the expressions for H(x) and G(x) are 

H(x) = \[cPP^-c^+vP'^y (9) 

{X> 4 F 2 2 FQ 4 F 2 2 Q 4 F 2 2 + 4 F 2 ' 1 j 

the primes meaning derivatives with respect to x. These equations were studied in |]J 
only under the simplifying assumption ^ = g(x) + h(y), which obviously corresponds to 
the vanishing of the function H. Given that this case was completed in the mentioned 
paper, we can assume from now on that H does not vanish identically. In the next 
section we are going to find the general solution of the system of partial differential 
equations (||) for the function *&(x, y) when H ^ 0. 

3 Solution of the system of equations for ^. 

In order to solve the system (^), let us find all the third-order partial derivatives of ^ 
with respect to the variables x and y. Deriving the first of the equations with respect 
to x and with respect to y we immediately get, respectively 

d xxy V = H'V + Hd x *, d xyy m = Hd y V, (11) 

and considering the derivatives of the second equation we find, after making use of ([TT|) 

d xxx m = Hd y m + Gd x m + G'V, d yyy V = H'V + Hd x m - Gd y m. (12) 

We can now impose integrability conditions on this set of partial differential equations. 
For instance, the difference of the partial derivative of the first equation in ( |Tl"D with 



respect to x with the partial derivative of the first equation in (|12"D with respect to y 
gives the following first order equation 

H"m + 2H'd x m = G'd y ty, (13) 

where the original system @ has been taken into account. It is easy to see that the 
integrability condition for the other two equations is given again by (|I3]) . 

If the function G is a constant, then ( [I3|) implies that either H is also a constant 
or ^ takes the form of a product of a function of x and a function of y. From ( |T3"D it 
also follows that when if is a constant then G(x) must also be constant (because the 
other possibility d y ^/ = needs not be considered as discussed above). Thus, let us 
first study the case when both H and G are constants. Since we can assume H ^ 0, 
the original system (H) implies the following homogeneous equation 

G 

<KA< - <), l; /\> - — d xy ^ = 
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which can be immediately solved (given that G and H are constants) to give 



a 2 H 



x 



where fx and f 2 are arbitrary functions of their arguments and the positive constant a 
is defined by 



-G + VG 2 + 4ff 2 
2H 2 



With this form for the function the second equation in 
into the following two equations 



l - (G + VG 2 + AH 2 ) . (14) 
can be trivially separated 



h = -^ [fi + 



a \ 
2HJ 



a = /*v(a-£ 



(15) 



where a is a separation constant and dots indicate derivatives with respect to their 
respective arguments. The redefinitions fi + ^ — ► fi and f 2 — ^ — ► / 2 , which leave 
invariant its sum & = fi + f%, allow us to set the constant a = 0. Thus, the solution of 
the system ([15]) is clearly 

Solution 1: * = Ae aHy+ « + Be~( aHy+ ^ + C sin ^ - atfa; + ar ) 

where A, B, C are arbitrary constants and Xq IS 8j superfluous constant that can be 
reabsorbed into the coordinates. Let us remember that G is obtained in terms of H 
and a from Q141) . 

We shall momentarily discard the case when G is a constant and ^ is a product of a 
function of a; and a function of y because it is contained in a more general situation we 
will meet later. Let us then consider the general case when G' ^ (and then H' ^ 0). 
Deriving equation fllTf ) with respect to x and using the first equation in (|j) we get 



IE' 



'2ir 

G' 



+ 



H" 
G 7 



* = 0. 



(16) 



The nice property of this partial differential equation is that all the coefficients and all 
the derivatives involve only the variable x. Thus, its general solution is 



* = +A 2 (y)v 2 (x) 



(17) 



where A% and A 2 are arbitrary integration coefficients which can depend on y, and V\ and 
t>2 are two linearly independent solutions (so that in particular they are non-vanishing) 
of the ordinary second order differential equation given by the same expression (flB] ) 
with ordinary derivatives substituted for the partial derivatives. Our aim now is to 
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find all the solutions of the original system ([8]) with the form for the function \1/ given 
in ([171) . It is obvious that we will also find all the solutions which are products of a 
function of y and a function of x (say by imposing A 2 (y) = in flTT|)), and thus the 
case G' = with H' ^ we had left aside is also contained in this study. 

The first equation in the system (|8]) for a function ^ of the type ([H]) takes the form 

Aiv[ + A 2 v' 2 - A x Ev x - A 2 Hv 2 = 0, (18) 

where dots mean derivative with respect to y from now on. The lefthand side of this 
equation is a linear combination of the four functions of y, Ai, A 2 , Ai and A 2 , with 
coefficients not depending on y. Thus, we can study this equation by considering 
all the possibilities of linear independence of these four functions. First of all, it is 
obvious that the case in which all four functions of y are linearly independent need 
not be considered as it implies the vanishing of each coefficient in (|H|) which is clearly 
impossible. Thus, we can move to the case in which Ai and A 2 are linearly independent 
and rank-f^!, A 2 , Ai, A 2 } = 3, so that without loss of generality we can consider that 
A 2 is linearly independent of A x and A 2 and write Ai = f3\A 2 + f3 2 Ax + P 3 A 2 , where f3x, 
P 2 , p 3 are constants. The equation (|T8| ) becomes 

A x (/3 2 v[ - Hvx) + A 2 (f3 3 v[ - Hv 2 ) + A 2 (/3 lV [ + v' 2 ) = 0, 

so that each term between brackets must vanish. From the coefficient multiplying A\ 
and H ^ we have that (3 2 ^ and thus v[ = Hvi/f3 2 . Substituting this value 
of v [ into the second term in brackets and dropping a global factor H we obviously 
get v 2 = P3V1/ j3 2 so that the two functions v\ and v 2 are linearly dependent which is 
impossible. Thus, this possibility gives no solution for the function \1/ either. 

Let us now consider the next possibility given by A\ and A 2 linearly independent 
and rank{^4 1; A 2 , Ai, A 2 } = 2. Thus, each function A\ and A 2 is a linear combination 
with constant coefficients of Ai and A 2 : 

A, = p\A x + 0[A 2 , A 2 = $\A X + (3 2 2 A 2 . 

We can redefine the two functions A\ and A 2 using arbitrary non-singular linear com- 
binations of themselves because this does not change the form ( |ITD of \1/ (obviously, we 
must redefine the two functions v\ and v 2 appropriately). Under such a change, the 
coefficients f3 h a (a, b = 1,2) transform like an endomorphism that can be rewritten in its 
canonical form. Consequently, the only two inequivalent possibilities that can arise are 

Ai = /M! \ Ai = (3 A x 

A 2 = (3 2 A 2 . J ' A 2 = Ax + (3A 2 

In the second case, equation ( |T8"D takes the form 

Ax (Pv[ - Hvx + v' 2 ) + A 2 ((3v' 2 - Hv 2 ) = 




which implies that each term enclosed in round brackets must vanish. Again (3 must 
be non- vanishing so that we can set (3 = -. The solutions for v\ and v 2 are trivially 

Vl = r(x)e ar{x \ v 2 = -\e ar{x \ with r' = H, 

a 2 

while the solutions for A\ and A 2 are 

A 1 = Ae%, A 2 = A(y - y )e*, 



where A ^ and dq are constants. Obviously, we can reabsorbe y$ into the variable y. 
Finally by imposing the second equation in (||) and using the above results we obtain 
the following two equations 



r" + 2ar' 2 + ^ = 0, G + a 2 r' 2 + A = . (20) 
cr a 2 

The first equation gives r(x) explicitly (and thereby the function H(x) after derivation), 
and the second gives G(x). Thus, both of them are differential relations for the unknown 
metric coefficients F(x), Q(x) and P(x). Summing up, we have 

Solution 2: * = Ae ar{x)+ « (r(x) - with H = r' 

and the two relations (|20|) must hold. 

We must now consider the possibility given by the first set of differential relations 
in fl!9|) . Using these equations, the relation ( jI8|) takes the form 

A, ((3 lV [ - Hv x ) + A 2 {(3 2 v' 2 - Hv 2 ) = 0, 

so that again each term in brackets must vanish. As neither V\ nor v 2 can be zero, 
we have that the constants f3\ and f3 2 are non-vanishing so that we can rename them 
respectively as ^ and The equations for A\ A 2 , v[ and v' 2 allow us to write the 
final form of \1/ immediately as 



Solution 3: * = Ae air{x)+ ^ + Be a2r{x)+ % , H = r' 



(with oi\ ^ a 2 in order to have A\ linearly independent of A 2 ). The second equation of 
([8]) is equivalent now to the two equations: 

n \ < i x /2 , «i + «2 n n '2,1 1 1 n / 01 n 

r +(01 + 0:2)^ H ^ — 5 — = 0, G + a\_a 2 r H — 5- H — ? H = 0. (21) 

afal a i ^2 «i«2 

It is convenient to note here that the constants ot\, a 2 , A and B can be complex still 
giving a real function \E'. This can be accomplished only when ot\ = a + i/3, a 2 = a%, 
A = B where the bar means complex conjugation. In this case the function \l/ becomes 

m = Ae ar(x)+ ^ y sin ( (3r(x) - f v + x ) 

\ a 2 + p 2 J 



9 



where xq is a constant that can be made zero by redefining the variable x. 

Let us finally consider the case when the two functions A\(y) and A 2 {y) are linearly 
dependent. Due to the symmetry in the indexes 1 and 2 we can choose A 2 proportional 
to Ax, A 2 (y) = XAi(y) where A is an arbitrary constant. The form of the function \1/ is 
now 

*(z,y) = A x {y) [vt(x) + Xv 2 (x)} = A x (y)X{x) 

so that it is the product of a function of only x and a function of only y. The first 
equation in (H) is now 

X' A x = H(x)XA l =► Ai = -A u X' = aHX 

a 

for some non- vanishing separation constant a. The second equation of the system (H) 
becomes (after dropping a global factor 

G(x) - ar" - aV 2 + ^ = 0. (22) 

a 1 

Thus, the final possibility for the system (^) is given by 

Solution 4: m = Ae ar{x)+ « , r' = H 



together with the differential equation (|22|). 



4 Analysis of the remaining Einstein equation 

Having found all the possible solutions of the system (Q), we must now impose the 
quadratic equation ([5]) on each one of these possibilities in order to find which of them 
give rise to a perfect-fluid solution of Einstein's equations. In order to deal with the 
mentioned quadratic equation it is worth considering it as the vanishing of the norm of 
a three-dimensional vector with components 

S = (Sqo + S33, Sn — S22, Soi) (23) 

in a three-dimensional vector space endowed with a scalar product given by 

/ i 
\ 

V -1 
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This quadratic form has signature (1,-1,-1), and therefore, equation (|) states that 
the vector (|23| ) is null in this indefinite scalar product. Now, evaluating the components 
of the Einstein tensor of the metric (HT) we find 



Son + S; 



33 



+ R 2 (x)V, x + m m + 2c*,j, 



Sn - S33 = * Mi(x)* + M 2 {x)^ )X - 2tf . 



01 



where the expressions for Fi, F2, Mi, M2, Li, L2 and L3 are 





F /2 
~ p~2~ + 


F'Q' 
FQ 




R 2 = - 


F' 

f~' 




2 P 2 Q 2 

° F 2 


+ c 2 


2 Q 2 




F' 

_ f"' 


F' 1 
F 2 


3 P' Q ') 
2 Q 




,Q 

L 2 = P-, 
F' 


L 3 = - 


-f(cP + v 



Let us now study the consequences of equation @ when the function \l/ takes each one 
of its possible forms. Let us begin by what we have named Solution 1. 

Solution 1 As convenient notation let us define the following three functions of the 
two variables x and y 



k\(x, y) = e aHy+ a , k 2 (x, y) = sin y aHxj, k 3 (x, y) = cos y aHxj, 

so that \I/ is written as \I/ = Afci + B/k\ + Ck 2 . If in this expression C vanished, the 
solution would in fact belong to Solution 3. Similarly, if both A and B vanished, the 
solution would be also contained in Solution 3. In consequence, we must have C ^ 
and at least one of the A and B non-zero. By using the change a — > — a we can obviously 
assume A and C non-vanishing in what follows without loss of generality . 
By taking the derivatives of the function the vector ( |2"3"D becomes 



^ = AktXt + ^-X 2 + Ck 2 X 3 + Ck 3 X 4 , 
w k\ 

where the four vectors Xj = 1, 2, 3, 4) depend only on x and are given by 

X\ = ( Ri + — + 2a 2 H 2 + 2caH, M 1 + — - 2a 2 H 2 , L x + — + aHL 3 ) . 
\ a a a J 



(24) 
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X 2 = B ( R x - — + 2a 2 H 2 - 2caH, M 1 - — - 2a 2 H 2 , L x - — - aHL 3 ) , 
\ a a a J 

X 3 = M l + ^,LX 

X A = ( — - aHR 2 , -aHM 2 , —aHL 2 + — 
\ a a 

and we have included B in X 2 because B is the only constant that can vanish. Taking 
into account that the only vanishing linear combination of products of two functions in 
the set {k\, 1/&1, k 2 , k 3 } is hx/ki — k\ — k 2 = 0, the vanishing of the norm of the vector 
(p4]) implies 

X\ ■ X\ = X± ■ X 3 = X\ ■ X4 = X 2 ■ X 2 = X 2 ■ X3 = X 2 ■ X4 = X3 ■ X4 = 0, 
X 3 ■ X 3 - X 4 ■ X 4 = 0, 2AX X ■ X 2 + C 2 X 4 ■ X 4 = 0. 

If X\ ■ X 2 7^ then the vector X4 would be non-null and jXi, X 2 , xA would be a basis 

of the vector space. From the first set of equations above it would follow that X 3 should 
be the zero vector. Then, the second equation would be incompatible. Thus, we must 
necessarily have X\ ■ X 2 = and every scalar product of each vector Xi with each Xj 
must vanish. Therefore, all these four vectors are proportional to some null vector and, 
as an obvious consequence, we have that the linear combination 



X\ — X 3 H — — — -X A 
a z H z 



a 2 H 2 + 4 



a 2 



aH' ' all) 



has to be also a null vector, that is to say 



A Ac L 3 2 

from which we learn that c cannot vanish. Consequently, Case III is impossible and we 
can set v — 0, so that from the definition of L 3 we have L 3 = —c^-. The equation above 
shows that L 3 must be constant so that every vector is proportional to the constant 
null vector (2 + 2c/aH, —2, L 3 /aH). In particular, from the proportionality of X 4 with 
this vector and the fact that H is actually a constant in Solution 1, it immediately 
follows that y- and ^ are also constants. Now, it is very easy to find the most general 
solution of Einstein's equations in this case, given by the following functions 

§- = -, P = e^r, — = lc =-2aH, * = Ae aHy+ ^ + C sin ( V - - aHx) . (25) 
Q a Q F \a J 

where e is a sign and 7 an arbitrary constant. The Einstein tensor for the metric (^) 
with these explicit forms for the functions F, Q, P and \l/ takes the following form 



5*00 + £22 — Sn — S 22 — —eSoi, S 22 — S 33 — 3C 2 ( a 2 H 2 H — - ) , S 23 — 0. 



a 
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Thus, it is clear that the Einstein field equations (|5|), (|6]) and (|7|) are satisfied. However 
this metric does not represent a perfect fluid because the fluid velocity vector is not de- 
fined. The reason is that the quadratic equation (|^) states exactly that the determinant 
of the 2x2 matrix 

I Sqq + S22 S01 \ 
\ $oi Su — S22 J 

vanishes. This is in turn equivalent to the existence of two functions Uq and U\ such 
that 

( S 00 + S22 Soi \ _ ( U$ U U X \ 



a 



2 



V $01 Sn — S22 J \ UqUi U 1 

where a is the sign of the component Soo+Sn (if it vanishes then a is the sign of — 5*22) 
if both are zero then the whole matrix is vanishing). When these two functions U$ and 
U\ satisfy 

ul -Uf>0 

we can normalise this "vector" and define the two components uo and u\ of the fluid 
velocity vector as 

aUl = (p + p) ul, aU^ = (p + p)uj, with Uq — u\ = 1. 

These expressions also define p+p. However, other possibilities may also happen. First, 
there exist solutions of the quadratic equation (H) which satisfy 

ul - ul < 0. 

This case would correspond to a perfect-fluid-like Einstein tensor but with the "fluid" 
velocity being spacelike and thus it has no physical interest as it does not describe any 
known matter. Finally, there exist solutions satisfying 

^0 ~~ Ui = 0, <^=^> Soo + S22 = Sn — S l 22- 
In this case the Einstein tensor is of the form 

S a p = K a Kp+ pg a0 

where K is a null vector lying in the two plane spanned by the Killings d t and d^. 
Thus, this solution represents a radiative fluid with pressure. So, we conclude that 
the quadratic equation (||) contains solutions which do not represent perfect-fluids and 
that this must be checked after having found the solution. It is clear now that the 
solution ( Po] ) belongs to this class of radiative fluids with pressure and thus will not be 
considered here any further. We can now move to the following solution given by 
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Solutions 2 and 3 Let us define the function 

y 

W(x,y) = r(x) - 

cr 

so that $ for Solution 2 is written as \I/ = Ae ar ^ + ^W. The calculation of the vector 
(f2~3|) gives now S = *ffAe ar ^ + * WY\ + Y 2 where the two vectors Yi(x) and Y 2 (x) are 
given by the expressions 

Y x = (r x + ar'R 2 + — + — , M x + ar'M 2 - — , L x + ar'L 2 + — 
V or a ai a 

f 2 = fr'iZa - ^ - \ r ' M 2 + 4, r'L 2 - ^§ 
V a 13 cr a 13 cr 

From the fact that W depends on y while the two vectors Y\ and Y 2 involve only 
functions of x it follows that the vanishing of the norm of the vector S is equivalent to 
the three equations: 

Y x .£ = 0, n-y 2 = 0, F 2 -F 2 = 0. (26) 

This is a system of ordinary differential equations to be satisfied by the unknown metric 
coefficients. At this stage, however, it is convenient to study the possibility given by 
Solution 3 rather than studying this system of equations because, as we will see, the 
two systems of equations for Solutions 2 and 3 can be treated together. 
For Solution 3 it is worth defining the two functions 



y 1 (x,y) = e airm+ ^, ^2(x,y) = e 



i 



so that \I> is written as \I/ = A^i + B^ 2 , where A and B are non- vanishing constants 
and «i 7^ a 2 . Taking the partial derivatives of the function the vector (|23"D takes the 



form S = ^ 



A^\Z\ + B^ 2 Z 2 , where the two vectors Z±(x) and Z 2 [x) are 



2 2c 2 L 

Z\ = [Ri + a 1 r'R 2 + + — , M 1 + a 1 r'M 2 5 , L x + a x r'L 2 + — 

V af cti af a>i 



, 2c ,2 L 3 

Z" 2 = LRj + a 2 r'i2 2 + + — , Mi + a 2 r'M 2 5 , Li + a 2 r'L 2 + — 

\ a 2 a 2 a 2 a 2 



As , ^1^2 an d ^ 2 are linearly independent functions of y, it follows that the quadratic 
equation ([5]) is equivalent to the vanishing of all the scalar products of Z\ and Z 2 so that 
both vectors must be proportional to some null vector. Of course any linear combination 
of them has also the same property. It is convenient to consider the linear combination 
given by 

Zi - Z 2 = «i - a 2 r i? 2 5-5 , r M 2 H 5-= — , r L 2 . 

y afa 2 a x a 2 afa 2 a x a 2 I 
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As ot\ — a 2 is non-vanishing, we can drop this global factor and the full set of Einstein's 
field equations to be satisfied in this case are given by 

r" + (<*! + a 2 ) r' 2 + = 0, (27) 

afar; 

H = r', (28) 

G + a x a 2 r' 2 + (ai t? )2 - — = 0, (29) 
afa 2 a\a 2 

(r x + a ir 'R 2 + 4t + — ) (m x + a x r'M 2 - 4^ - (h x + a x r'L 2 + — ) = 0, (30) 
\ a( aij \ afj \ a x J 

( J* + a^ifc + A + -1 + 2^+^) + (r'ifc - 2^i±^- 

2c \ / . , , , , 2\ _ / _ \ / , , L 



Mi + ai r'M 2 - _ ] - 2 f Li + «ir'L 2 + — ) fr'L 2 = 0, (31) 



ol\OLi) \ a{) \ ol\) V aia2/ 

U, - 2^ - -^-) f/M 2 4- 2*+*) - (/L 2 - ^) = 0. (32) 
At the sight of this system of ordinary differential equations, we note that although 



Solution 3 is characterized by the property a\ ^ a 2 , the system (|27|)-(|32"D itself can 
be particularized to a\ = a 2 without problem (this is due to the fact that we have 
dropped a global factor a± — a 2 above) and gives exactly the system (|26D which had to 
be satisfied for the previous Solution 2. Thus, by studying the compatibility of system 
(p7| )- (j32]) we are considering both Solutions 2 and 3 together and thereby we have 
reduced considerably the amount of work to be done. By assuming that the resulting 
spacetime has a well-defined axis of symmetry we can expand the functions involved 
in this system around the axis of symmetry and check the compatibility of this system 
of differential equations. In the Appendix we show that that the system (|27|)-(|3"2"D has 
no solution giving rise to a spacetime with regular axis of symmetry. Let us emphasize 
that this assumption of regular axis is not an additional condition in our work because 
the restriction of the possible Lie algebras to the four cases listed in the Introduction is 
based on a geometrical result proven in which makes essential use of this regularity 
property. 

Solution 4 This last case is the most interesting one because it gives rise to a family 
of stationary non-convective differentially rotating perfect-fluid solutions. In fact, from 
the explicit form of \1/ = Ae ar ^ + ^ the vector ( p3f) takes the form 

S = A^ 2 (r x + ar'R 2 + + — , M x + ar'M 2 - L x + ar'L 2 + —) , 

so that the condition that this vector is null gives only one condition for the three 
functions F, Q and P. Besides, we know that for this Solution 4, H and G had to 
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satisfy only the two relations H = r' and (|22|), with no relation involving only r. Thus, 
the Einstein field equations reduce to a system of three ordinary differential equations 
for the four unknowns r, F, Q and P. We devote the following section to consider these 
perfect-fluid solutions. 



5 The solutions 

In order to write down explicitly the field equations and the density and pressure for 
the family just found, let us rewrite the line-element in a separate way for the Cases II 
and III. It is obvious that the metric for the Case II can be rewritten as 



ds 2 



D 2bcy 



F(x) e~ cy dt + P{x)d<j> 



Fix) 



1 7 2 1 

ax + 



L 2 (x 



L 2 (x 



■dy 2 



(33) 



where b is a constant and the functions F, Q and P are different from those used in the 
previous sections. We do not use different symbols for this metric coefficients in order 
to avoid using too much notation. In terms of these functions, the field equations to be 
solved read 



F' F 2 T 



IF* 
2 F 2 
1 / F' 



F Q + L Q 2Q 2 



0. 

c 2 p2 F 2 



^ + 

Q 2 Q 2 



+ c 2 (2b 2 



IF' Q' 



V 



V Q' 1 F 2 



, P 2 F 2 

~Q 2 



- — + - — ^ + - + -^r + -^P' Z + 2bc 2 (b-l) + c 



2 \ F 



2 F Q 



L 



LQ 2Q 2 



2\Fl 9FO \l) LO + ?0 2 O °\ 



2 F Q \L 



L Q 2Q 
F 2 



Q 



(34) 
(35) 



Q 2 



Hp' + hp" _ \9Lp> + C 2 P t b _ r 

F 2 2Q y 



1 2 



• (36) 



The explicit expressions for the density and pressure are simply 



P 



e 2bcy 
P = - 



IF'' 



1 F' Q' V Q' 1 F 2 . 2 



4^ + C 



AF 2+ 2F Q L Q + AQ^ 



L 2 



,2bcy 



'IT 



L Q 2 Q Q 2 



+ c 2 2b 



while the fluid velocity vector of this solution reads 



u 



~ bcy A( 



x) 



d_ 
dt 



+ R(x) 



d 
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where A(x) and R(x) are rather long expressions depending on the functions F, Q, P 
and L. We note that the solution is differentially rotating except when R(x) = 0, which 
corresponds to the case when both the second and third terms in square brackets in the 
last equation are vanishing. 

The system of equations (|34|)-(|36|) is formed by three ordinary differential equations 
for four unknowns so that, in principle, it contains an arbitrary function in its general 
solution. However, the third equation in the system is quadratic so that the compatibil- 
ity of the system is not obvious at first sight. In order to see explicitly that the system 
is indeed compatible and contains an arbitrary function, let us define a new function 
M(x) by 

p 

Thus, we are assuming P / (the case P = will be dealt with below). Notice that 
when P is non-zero we can always make it positive (locally) by using the coordinate 
change <fi — > —(f). Now, we define three new functions u(x), s(x) and w(x) by 

Q' P' V 

We will consider M as the arbitrary function so that the system gives differential equa- 
tions for u, s and w. First of all we notice that the case M = 1 gives a solution which is 
not a perfect fluid (see the discussion above) and, therefore, we can assume M^l from 
now on. From equation ( ffif ) we obtain u and substitute it into the other two equations. 
Then, fl35|) allows us to obtain s' in terms of the rest of the functions. Substituting into 
(|36|) we obtain an expression of the form 

w' 2 (M - 1) (46 2 - 46 + M) + u/Ei (M, M'M", w, s) + S 2 (M, M', M", w, s) = 0. 

When M = 46 — 4b 2 the system (p5|)-(p6|) is linear in their two highest derivatives (w', s') 
and therefore the existence theorems apply. In the more general case M ^ 46 — 46 2 , 
we have to solve w' from this quadratic equation to obtain a system which is linear 
in their highest derivatives. This can be done unless the discriminant of this equation 
is negative. Evaluating this discriminant we notice that contains both positive and 
negative terms so that it is not negative definite for all M and all initial conditions. In 
general, there will be functions M and initial conditions for n and s (as well as for the 
constant 6) which make the discriminant positive thus giving solutions to the system. 
Just to give an explicit very particular example (they are trivial to construct), the value 
6=1, the initial conditions w \ x =o= 0, s \ x= q and any function M with M \ x =o= 1/2, 
M' I x =o= and M" | x =o= give a positive discriminant and produce a well-defined 
solution. Thus, we conclude that, even though we cannot prescribe the function M 
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freely, the general solution of the system (|34T)-(|3"6|) contains an arbitrary function M 
belonging to a suitable open subset of the whole class of C 2 functions. 

Let us now consider what happens when P = (which corresponds to the static 
limit of the family). It is indeed very simple to show that the equations can be rewritten 
in this case in the form 



u 

— = w{x) + 0, 



— = 2 26-1 _ = _ + — _ — — w 

F Q w b — 1 



w 
w 



86 2 -46-1 , (26 2 
V + 



6-1 



(6-iy 



(26-l) 2 u> 2 + (6-l)"c 



2 J2 



(Mi) 

(37) 



(the case w = is slightly different and can be integrated completely to give a solution 
which contains more symmetry). The energy-momentum contents for the solutions of 
(|37j) is that of a fluid with velocity vector proportional to d t and density and pressure 
given by 



P 



L 2 (6 + 1) 
(6-1) 



D 2bcy 



l)w' + (2b-lfw 2 + (b-l) 2 c 



\2 2 



P 



b + r 



so that it satisfies a barotropic equation of state. 

We consider now the Case III. The metric takes the form 



ds 2 = e 2bvy 



-F(x)dt 2 + 



F{x) 



+ (P(x) - vy)dt 



L 2 (x) 



dx 2 + dy' 



(38) 



where the functions F and Q have been again redefined by introducing L(x). Einstein's 
equations with these new functions are 



L AF 2 



0. 



1F ,Z F'Q' L'Q' 1Q 2 , 2 Q" , 2 2 v Q 2 

— + 2 — — P' 4- — I- 26 v A — = 

' FQ + LQ 2F 2 + Q + ° + 2F 2 U ' 



2 F 2 

F 
'V_ 



(39) 
(40) 



I F'Q 



V Q' 1 Q 2 



- - + + -r + ~ - ~P 



1/2 



2 F Q 

y_q__^ 

LQ F 2 



L 



L Q IF 2 ' 



26V 



1 F" 1 F' Q' 
YY ~ 2~F~Q 



P 



2 F 2 



F 2 



F 2 2 Q 



(41) 
(42) 



The fluid velocity vector of this solution reads 



u 



X) 



dt 



+ [vy + R{x)\ 
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where, as before, A(x) and R(x) are expressions depending on the four functions F, Q, 
P and L. We note that in this case the solution is always differentially rotating (its 
rigidly rotating limit must satisfy v — and thus it actually belongs to the case when 
the conformal Lie algebra is abelian). The expressions for the density and pressure for 
this solution read 



V 



D 2bvy 



1 F' 1 1 F' Q' 



F® itf 2 v^tf 
L Q 4F 2 4 F 2 



D 2bvy 




,F_q_ 

L Q 



191 

2 Q 



+ 2—^- + + 2b 2 v 2 + v 2 



F 2 



The same considerations we made regarding the existence of solutions for the qua- 
dratic set of differential equations (|34"l)-(j36l) hold in this case. Given that the equations 
are simpler, the analysis is easier than before. We simply need to define a new function 

N(x) by 

Q = FN 

and regard N as the arbitrary function in the system. Introducing two new functions 

s(x) and w(x) by 



F' 
~F 



F 



w. 



we notice that (^9|) allows us to obtain P' in terms of and w. Substituting into 
the other two equations we find that ( f4"0D contains s' linearly and no other derivatives 
(except for N). Substituting s' obtained from this equation into the last one we find, 
similarly clS cLS before, 



46 2 



( 1 + W2 j w' z + Ei (N, N', N" , w, s) w' + E 2 (JV, N', N", w, s) = 0. 

Evaluating the discriminant of this quadratic equation it can be seen that it is not 
negative definite. Thus, there exist initial conditions and functions iV which make it 
positive. A trivial example is given now by s | x =o= 0, w | x =o= 0, 6 = 1/2 and any 
function satisfying A^ \ x =o— 1> N' \ x =o— and A^" \ x =o= 0. Therefore, the general 
solution of the system (|39D-(f4"l~D will contain an arbitrary function belonging to some 
open non-empty subset of the class of C 2 functions. 

It is interesting to note also that the conformal Killing of all these solutions (both 
for algebra II and algebra III) is in fact homothetic as it satisfies 



£%9w = 2bcg 



(Case II), 



^k9^ = 2 bvg 



(Case III), 
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so that the scale factor is a constant. Thus, the barotropic equation of state of the 
perfect fluid, in case it exists, must be a linear relation of the form p = (7 + 1) p for 
some constant 7 (see J4|). As the families of solutions contain an arbitrary function, 
there always exist solutions satisfying this barotropic equation of state. 

Let us finally summarize the main results of this paper and Ref. JI| in the following 
table where the solutions of Einstein's field equations for stationary and axisymmetric 
space-times with a (non-isometric) conformal motion (proper for the Abelian Case and 
Case I, but homothetic for the remaining Cases II and III) and filled with a non- 
convective perfect fluid are written down. 



Abelian Case 




Schwarzschild interior solution 
Solution in || 


Case I 


Wo 
** = o 


Schwarzschild interior solution 
Solution in M 


Case II 


\P j, = -fee* 


Solution (j33|) 


Case III 




Solution ([§) 



Corrigendum to paper 

Let us end this paper by adding a corrigendum to our previous paper The general 
solution for the Case I was explicitly found there. However, the expression for the fluid 
velocity vector (the expression following formula (27) in M) is wrong and should read 

1 d 

(3^H 2 - H 2 u 2 X 2 dt 

which implies that the rotation Q of the fluid is constant. Thus the general perfect 
fluid solution in Case I is rigidly rotating contrarily to what was stated in [|TJ. As a 
consequence, we have that this solution must be contained in the general Wahlquist 
family || which is the general solution (together with its limit cases [000) for Pctrov 
type D stationary and axisymmetric rigidly rotating perfect fluids with equation of 
state p + 3p =const. (see 0). 
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6 Appendix 

The aim of this Appendix is to show that the overdetermined system of ordinary differ- 
ential equations fl2"7|)-(|3"2]) has no solution representing an axially symmetric spacetime 
with a well-defined axis of symmetry. First of all, we need to find the conditions on the 
metric functions F, P and Q such that the metric (|j) has well-defined axial symme- 
try. Thus, we must impose the existence of a two-dimensional surface where the axial 
Killing vector vanishes and where the so-called elementary flatness condition (see 
e.g. |10[) is fulfilled. A trivial calculation shows that these conditions are equivalent to 
the existence of a finite value of x (which can be redefined to x = using a coordinate 
change x — ► x + const.) where 

Q \ x =o= 0, F \ x=0 = 1 and Q' \ x=0 = 1 (43) 

(in order to find these expressions we have also made use of the coordinate freedom of 
rescaling the coordinate t). We must now expand the metric functions around x = 
and show that the system of ordinary equations (f2"T|)-(|3"2|) is incompatible. The first 
thing we notice is that the equation fl27|) for r can be completely integrated to give 

1 

r = tan 

(X\Ol2 



X — Xq J 



if oti + a 2 7^ 0, 
r' = const ^ if Oi\ + a 2 = 0. 

In the second case we assume the constant r' to be non-zero because otherwise the 
equation ( p8|) gives H = which belongs to another case. Now, a very simple analysis 
of the function \I/ (which is given by Solution 2 or Solution 3) shows that r' must be 
finite at x = in order to obtain a non-singular metric at the axis of symmetry x = 0. 
This implies that a ^° 2 x ^ ± (| + nnj , (n being an integer). This result is necessary 
in order to consider properly the different powers in x of the equations. 

We must now consider separately the case III from the case II. The first one (c = 
0,v 7^ 0) is simpler to study and the proof is as follows. First, we obtain 

_ 4r' F 2 
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from the equation fl28|) and substitute this expression everywhere so that only the 
functions F, Q and r' are present in the equations. We define a new function w(x) by 
Q = xw(x) so that w(0) = 1 (see (|43|)). The coefficient in x~ 2 of equation fl32|) gives 
simply 



16 ,4 , 

— r | a . =0 = 0, 



v 



which implies from the expression for r' above that Xq = and ati + a 2 ^ 0. Now, the 
coefficients in a; -1 of equation and (p?TD read, respectively 



-- [F' U =0 -2m/ | x=0 ] = 0, t F ' U=o -4^' U=o] = 

2 aran 



so that obviously F' \ x =a= and w' \ x =o= 0. Using these expressions the coefficient in 
x~ 2 of (Bol) reads 



16 («i + «2 



9 4 4 



) 2 

- = 



which is impossible and we have shown that no regular solutions around the axis of 
symmetry exist for the case III. 

Regarding the case II (c ^ 0, v — 0), the proof is more involved because the equations 
are more complicated. The idea is, however, the same as in the previous case. Now it 
is convenient to define a new function M(x) by 

" " Q (44) 



,'M 

which allows us to write (]28 ) in the form 



P M-l \ 2 M Q c J V; 

(the subcase M = 1 can be easily shown to be impossible) and substitute ^ everywhere 
so that only the functions Q, M and r' appear in the system. We have to take into 
account, however, that we still do not know the behaviour of the function M near 
x = because the regularity on the axis of symmetry does not tell us anything about 
the function P (and therefore about M through its definition), except that QP must 
be finite at x = (because the norm of the timelike Killing d t must be finite and 
negative on the axis of symmetry) and, therefore, the function M must be regular at 
x = 0. Our next aim is to show that necessarily r' \ x= o= (or, equivalently, Xq = and 
ax + a 2 7^ 0). Let us suppose, on the contrary, that r' \ x =ot^ (but finite as discussed 
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above). Again we define the function w(x) by Q = xw(x) and substitute everywhere. 
The coefficient in x~ 2 of (|3T 



gives 



-M 



=o r 



,/2 



1=0" 







M 



x=0~ 







and using this information, the coefficient in x~ 3 of fl3~0|) gives M' | x= o= and then the 
coefficient in x~ 2 in (|3~TD gives M" | x= o= 0. Thus, we have that M has a zero of order 
n > 3 at x = which implies from its definition that P has a zero of order | — 1, 
but then the expression (|45D is impossible in the limit x = and we have proven that 
r'U=o=0. 

It is convenient to perform the change of variables 



u 



tan 



(c*i + a 2 ) 
a\a 2 



x 



so thar r' = -^^u and the axis of symmetry is located at u = 0. Now, expanding 
M and w around w = it is not difficult to show that the system of equations is 
incompatible. The calculations are, however, long and the use of computer algebra is 
needed. Moreover, a number of different subcases arise so that the proof is lengthy and 
cannot be reproduced here. Let us simply give the main indications for the proof. It is 
enough to restrict the analysis to the three simplest equations ( 
it is convenient to redefine the constants to 



, and (H) and 



«i = 



K-Ko 



a 2 = 



1 



and define a function n(x) through 

Q " 



(«i + a 2 ) n 



Finally, the system of equations takes a much simpler form by performing the following 
combination 



Equationl = 4r /2 (f£| - 
Equation2 = 2r /2 (|3l|) + r' 



Equation3 = (|32| ) 



dx 



r (cti - a 2 ) H 5-^5— (n 

afa 2 



The first of these equations contains M linearly and without derivatives of M. This 
allows us to substitute M from that equation into the two others so that we have two 
equations for the single function n(u). The analysis now is lengthy but trivial. 
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